Antiferromagnetism in two-dimensional quasicrystals 
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Quasiperiodic structures possess long range positional order, but are freed of constraints imposed 
by translational invariance. For spins interacting via Heisenberg couplings, one may expect therefore 
to find novel magnetic configurations in such structures. We have studied magnetic properties 
for simple two dimensional models, as a first step towards understanding experimentally studied 
magnetic quasicrystals such as the Zn-Mg-R(rare earth) compounds. We analyse properties of 
the antiferromagnetic ground state and magnon excitation modes for bipartite tilings such as the 
octagonal and Penrose tilings. In the absence of frustration, one has an inhomogeneous Neel-ordered 
ground state, with local quantum fluctuations of the local staggered magnetic order parameter. We 
study the spin wave spectrum and wavefunctions of such antiferromagnets within the linear spin 
wave approximation. Some results for spin-spin correlations in these systems are discussed. 

PACS numbers: 71. 23. Ft, 75.10.Jm, 75.10.-b 



Recent experimental investigations of rare-earth based 
magnetic quasicrystals raise the question of what types 
of magnetic order are expected to exist in such struc- 
tures. The Zn-Mg-R series (R: rare-earth) are an exam- 
ple of localized spins in a quasiperiodic geometry, where 
evidence has been found [1] for magnetic correlations at 
low temperatures resulting from predominantly antifer- 
romagnetic interactions between spins. We therefore con- 
sider the issue of finding novel quantum magnetic states 
in this type of system at low temperatures. This ques- 
tion will be addressed here for the simplest theoretical 
situation, in which there is a spin on every vertex of a 
perfect quasiperiodic tiling. Study of the experimentally 
relevant effects of diluting the system and of disorder, 
which are nontrivial, are reserved for subsequent studies 
[2] . We describe the structure of the ground state and the 
spectrum of excitations for antiferromagnetically coupled 
Heisenberg S = 1/2 spins for this simplest case. At zero 
temperature, despite the fact that thermal fluctuations 
are frozen out, the spins do not order perfectly, due to 
quantum fluctuations. We consider the case of 5= |,the 
case of greatest theoretical interest, as quantum effects 
are expected to be strongest. 

We consider spins located on vertices of two- 
dimensional tilings: the octagonal tiling formed from 
squares and 45° rhombuses, and the Penrose tiling 
formed from thick (72°) and thin (36°) rhombuses. The 
eight-fold symmetric octagonal tiling is a somewhat sim- 
pler structure compared to the Penrose tiling (five-fold 
symmetric), in terms of the local environments present. 
The latter is more closely related to the three dimen- 
sional icosahedral tiling relevant in the description of a 
large number of quasicrystalline alloys. In both cases, at 
T=0, since there is no frustration in the model, one ex- 
pects a Neel ordered ground state, with equal and oppo- 
sitely directed sublattice magnetizations. We have stud- 
ied the antiferromagnetic ground state and excitations 
of the quasiperiodic Heisenberg models using linear spin 



wave (LSW) theory. When possible, the LSW results 
will be compared with Quantum Monte Carlo (QMC) 
data that have been obtained for our models. The cal- 
culations involve a numerical diagonalization of Hamil- 
tonians for periodic approximants of different sizes, with 
a finite size scaling analysis for an extrapolation to the 
infinite quasiperiodic tiling. The approximants of the oc- 
tagonal tiling were obtained using the method described 
in [3], while the approximants of the Penrose tiling were 
obtained by a generalization of the method described for 
the 36 site Taylor approximant in 0]. 

Model of the quasiperiodic antiferromagnet 

We consider the nearest-neighbor antiferromagnetic 
spin-^ Heisenberg model, H 



JT.{i,j) ^i-^j where 



are pairs of sites linked by an edge. The antiferromag- 
netic coupling J > is the the same for all bonds. This 
model is bipartite, in that sites belong to one of two sub- 
lattices A and B, with a spin on sublattice A coupled to 
spins lying on sublattice B and vice- versa. This property 
ensures that there is no frustration, i.e., if one considers 
classical spin variables, the ground state is one for which 
all bonds are "satisfied" - with all the A-sublattice spins 
pointing in one direction and B-sublattice spins pointing 
in the opposite direction. The site-dependent staggered 
magnetization is defined by nisi = {^i^t) = S — dm 
where takes the values +1 and —1 on the A and B sub- 
lattices respectively. The correction term Srrisi arises due 
to quantum fluctuations. The global average of nisi over 
all sites i is denoted mg. Assuming that quantum fluctua- 
tions are small, one can replace the exact Hamiltonian by 
a linearized form. This LSW Hamiltonian is then diag- 
onalized using a generalized Bogoliubov transformation 
(as described in [5]), and the complete set of eigenvalues 
and eigenvectors are obtained. The local magnetizations, 
ground state energy and magnon modes are thus numer- 
ically determined. 
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FIG. 1: (left) Density of states for Taylor approximant 
(N= 11556 sites) as a function of energy (expressed in units of 
J), (right) Scaling of the ground state energy per site Eq/N 
with system size. 
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FIG. 3: (left) Staggered correlation function versus distance 
from a fixed central z = 3 site calculated for Taylor ap- 
proximant. (right) The same staggered correlation function 
versus inverse distance, showing several parallel linear trends 
(each color represents a different z-family) 





FIG. 2: (left) Real space distribution of staggered magneti- 
zations on a Penrose approximant (colorscale: smallest (red) 
to largest (blue)), (right) Staggered magnetizations projected 
onto the second window in perpendicular space (same color 
coding) 



The magnon spectrum. Ground state energy 

The density of states is plotted as a function of the en- 
ergy (expressed in units of J) in Fig. (1) (left) for the case 
of a Penrose tiling approximant. The main features are: 
a continuous low energy part, a peak corresponding to de- 
generate "loop-like" states E = 3J, and high energy 
bands corresponding to states centered around high co- 
ordination number sites. The integrated density of states 
is well- fitted by a quadratic energy dependence, showing 
that there is a linear dispersion of the magnon modes 
in this region of the spectrum. The effective spin wave 
velocity on the Penrose tiling is c = 1.08 J, and in the 
octagonal tiling, our estimated value is Cocta ^ 1.3 J. For 
comparison , in the square lattice, Cgq = V2J ^ 1.41 J on 
the square lattice (for edge length a = 1 and S = 1/2). 
The Penrose tiling spin wave velocity is smaller than 
in the octagonal tiling, which is in turn smaller than 
the value on the square lattice, reflecting the hierar- 
chy in the areal density of sites. The ground state en- 
ergy per site is found by extrapolating finite size results 
as shown in Fig. (1) (right) (the fit to the A^~^/^ depen- 
dence is expected from LSW theory) and has the value 
Eo/N = -0.643(0) in units of J. This is close to the val- 
ues found in the octagonal tiling and the square lattice 
(two cases where analytical calculations are possible), de- 
spite spatially different ground states, primarily because 
the average number of nearest neighbors is the same in 
all three cases. 



Wavef unctions. Spatial distribution of the mag- 
netization at T=0 

Wavefunctions are extended and show multifractal 
scaling at low energies. One-dimensional states are found 
at the degenerate energy E = 3J, while the narrower 
energy bands at the top of the spectrum correspond to 
wavefunctions localized on high coordination sites. From 
these wavefunctions, one can compute the order param- 
eter in the ground state and study spatial variations. 
Fig. (2) (left) represents how local magnetizations vary in 
space on a portion of the Penrose tiling. The color of the 
circles around each vertex varies from red (small magne- 
tization) to blue (high magnetization). We find that the 
distribution of magnetizations can be explained qualita- 
tively in terms of a simple local cluster approximation - 
the Heisenberg star approximation (see Jagannathan et 
al in [s*] ) . The symmetries hidden in this complex ground 
state are best shown by representing it in a perpendic- 
ular space projection. Fig. (2) (right) shows the magneti- 
zations of subset of sites that project into one of the four 
pentagon windows in perpendicular space to illustrate 
this point. One can see the grouping of colors accord- 
ing to the environment. Note: in this representation, the 
closer two sites are, the farther the distance out to which 
their environments are identical. 

Correlation function. Structure factor. One can 

define, for a given spin in the tiling So, the staggered 
spin-spin correlation function C{f) = (coe^SoS^). As in 
periodic systems, this function turns out to decay with 
the inverse of distance r = |r|, but it also depends on 
the type of sites involved. This is shown in the left-hand 
side figure in Figs. (3), which shows all the spin-spin cor- 
relations of a fixed three-fold site placed at the origin in 
a Taylor approximant. The figure shows two groups 
of sites - more highly correlated ones (corresponding to 
the red upper branch) and less correlated ones (the black 
lower branch). The correlations are grouped in families, 
each of which decays as the inverse distance from the 
central site. This is seen more clearly in the right-hand 
figure, where C{r) is plotted against 1/r. Each color 
corresponds to a site of different coordination number 
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(for example, green points correspond to z = 4). The 
quasiperiodic character of the antiferromagnetic ground 
state is also seen in the static structure factor S{q}. Our 
computations show the apparition of magnetic peaks cor- 
responding to half-integer indices [5] as expected on gen- 
eral theoretical grounds. 

Conclusions. The studies show a rich structure of the 
magnetic ground state and correlations in 2D quasiperi- 
odic antiferromagnets. Comparison of the average stag- 
gered magnetization in the octagonal tiling and in the 
Penrose tiling shows that the former has larger quantum 
fluctuations compared to the latter. Spin-spin correla- 
tions have a rich structure, with some features of power 
law decay in common with 2D periodic systems. Interest- 
ing questions remain for future investigation, as to quan- 
tum coherence and entanglement in such systems, and 
the effects of perturbations, including phason disorder. 
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